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INTRODUCTION 


Kerma  (kinetic  energy  released  in  material)  is  defined  as  the  amount  of 
energy  per  unit  mass  transferred  by  indirectly  ionizing  radiation  to  charged 
particles  in  a  volume  of  material.  Kerma  is  equal  to  radiation-absorbed  dose 
under  conditions  of  charged  particle  equilibrium,  provided  that  energy  loss 
due  to  bremsstrahlung  is  negligible.  Kerma  factors  [expressed  in  units  of  rad 
per  unit  radiation  fluence  (rad  •  cm2)]  are  functions  of  the  radiation  type 
(neutrons  or  gamma  rays),  radiation  energy,  and  material  irradiated.  Kerma 
factors  are  used  in  neutron  gamma  dosimetry  for  a  reactor  to  calculate 
radiation  dose  when  the  fluence  is  known.  Ratios  of  kerma  factors  are  used 
to  calculate  radiation  dose  in  one  material,  such  as  tissue,  when  the  radiation 
dose  is  known  or  has  been  measured  in  another  material,  such  as  an  ion 
chamber  (1). 

Irradiation  by  neutrons  or  gamma  rays  of  a  spectrum  of  energies  requires  the 
use  of  spectrum-averaged  kerma  factors  (SAKF's)  to  perform  dosimetry 
calculations.  Detailed  energy  spectra  were  recently  determined  (2,3)  for 
several  configurations  of  the  AFRR1  TRIG  A  reactor.  These  spectra  were 
used  in  the  present  report  to  calculate  SAKF’s  and  ratios  of  SAKF's  for  the 
following  materials  of  interest  in  reactor  neutron  gamma  dosimetry  with 
paired  ion  chambers: 

(Kt)g»  (Kj)n 


(KT/Kc)g,  (KT/KMg)G 


(KT/KAr)N>  (kt/kC02)N> 

and  (Kt/Kte  gas)N 


METHODS 


Neutron  and  Gamma  Spectra 

Neutron  and  gamma  spectra  for  18  configurations  of  the  AFRRI  TRIGA 
reactor  have  been  reported  (2,3).  Table  1  summarizes  the  configurations  for 
which  spectra  are  available,  and  Figure  1  shows  two  sample  spectra.  The 
spectra  were  reported  in  the  format  of  Oak  Ridge  Data  Library  Collection 
DLC-31  (4),  which  consists  of  a  37-energy-group  format  for  neutron  spectra 
and  a  21-energy-group  format  for  gamma  spectra.  Angular  distributions  of 
the  neutrons  were  reported  in  terms  of  "front"  and  "back"  directions  of 
incidence,  but  only  the  total  ("front"  plus  "back")  spectra  were  used  in  the 
calculations  presented  here. 


Gamma  and  neutron  SAKF's  for 
tissue 

Gamma  SAKF  ratios  for  carbon 
and  magnesium  (ion  chamber  wall 
materials) 

Neutron  SAKF  ratios  for  argon, 
carbon  dioxide,  and  tissue- 
equivalent  (TE)  gas  (ion  chamber 
gases) 
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Table  l .  Reactor  Configurations  for  Neutron  and  Gamma  Spectra 


Room 

Distance 
to  Core 
(cm) 

Configuration 

Average  Ener 

Neutron 

gy  (MeV) 

Gamma 

Neutron 

Spectra  * 

1 

50 

Unshielded 

1.46 

1.16 

1-D 

100 

Unshielded 

1.49 

1.15 

1-D,  3-D,  MEAS 

200 

Unshielded 

1.32 

1.05 

1-D 

300 

Unshielded 

1.16 

0.91 

1-D 

400 

Unshielded 

1.11 

0.78 

1-D 

500 

Unshielded 

0.89 

0.67 

1-D 

2 

50 

Unshielded 

1.35 

1.18 

1-D 

100 

Unshielded 

1.32 

1.13 

1-D,  3-D,  MEAS 

200 

Unshielded 

1.12 

0.95 

1-D 

300 

Unshielded 

0.93 

0.76 

1-D 

NA 

NA 

Pneumatic  tubes 

0.88 

0.83 

1-D 

1 

100 

7.5  cm  H2O 

1.99 

1.04 

1-D 

100 

30  cm  H2O 

3.71 

0.95 

1-D,  3-D 

100 

5  cm  Pb 

1.58 

1.65 

1-D,  3-D 

100 

5  cm  Pb  and  exer- 

cise  wheel 

1.05 

1.36 

1-D,  3-D,  MEAS 

100 

15  cm  Pb 

0.96 

0.92 

1-D.  3-D 

100 

15  cm  Pb  and  cave 

0.68 

1.80 

1-D,  3-D,  MEAS 

100 

15  cm  Pb  and  evlin- 

drical  phantom 

0.55 

1.13 

1-D.  3-D.  MEAS 

1-D,  3-D,  and  MEAS  refer  to  the  availability  of  neutron  spectral  data 
derived  from  one-dimensional  ANISN  calculations,  three-dimensional  Monte 
Carlo  calculations,  or  activation  foil  measurements,  respectively  (2,3). 


15  cm  IEA0  WITH  PHANTOM 
15  cm  LEAD 


NEUTRON  ENERGY  ItV) 


Figure  1.  Neutron  spectra  at  100  cm  from  the  AFRRI  TRIGA  reactor,  shielded  by  15-em  lead  wall. 
Solid  curve  shows  spectrum  free  in  air;  dashed  curve  shows  spectrum  at  midline  in  an  H-em- 
diameter  cylindrical  tissue-equivalent  phantom. 


The  focus  of  the  spectrum  calculations  of  references  2  and  3  was  on  neutron 
spectra;  thus  more  reliance  can  be  placed  on  these  than  on  the  gamma 
spectra.  In  all  cases,  one  dimensional  (1-D)  ANISN  calculations  of  the 
spectra  were  performed.  The  gamma  spectra  excluded  fission  and  fission- 
product  gamma  radiation  since  secondary  gamma  rays  were  presumed  to 
dominate  the  spectrum  (3).  Therefore,  the  shapes  of  the  gamma  spectra 
(which  depend  on  the  materials  irradiated)  are  probably  realistic,  but  the 
absolute  magnitudes  (which  depend  on  the  ealculational  model)  must  be 
viewed  with  caution.  For  the  neutron  spectra,  the  calculations  were  refined 
with  three-dimensional  (3-D)  MORSE  calculations  in  nine  cases,  and  verified 
by  activation  foil  measurements  in  five  of  the  more  complex  cases.  Thus  the 
neutron  spectra  represent  the  best  state-of-the-art  determinations  now 
available. 

To  indicate  the  degree  of  moderation  of  each  neutron  and  gamma  spectrum, 
average  spectrum  energies  were  computed  according  to  the  following  formula 
(and  are  listed  in  Table  1): 


E  4ME)  dE 
0.4  eV 

Eavg  ~  » 

j  «l’(  E)  dE 
0.4  eV 


where  E  and  Ej  represent  neutron  or  gamma  ray  energy,  and  <ME)  and  Oj 
represent  the  number  of  neutrons  or  gamma  rays  with  that  energy.  Note  that 
this  definition  of  average  energy  excludes  thermal  neutrons. 

In  addition  to  the  spectra  listed  in  Table  1,  SAKF's  were  also  calculated  for 
the  reference  spectra  given  in  DLC-31  (4).  These  spectra  were  prompt 
gamma-ray  fission  source,  neutron  spectra  for  fission  and  thermonuclear 
sources,  and  14-MeV  neutrons.  A  final  reference  spectrum  was  the  ENDF 
uranium-235  fission  spectrum,  which  was  used  as  the  source  term  for  the 
AFRRI  neutron  spectrum  calculations  (2).  While  none  of  these  spectra  apply 
directly  to  the  AFRRI  reactor  configurations,  they  were  included  in  the 
present  calculations  for  comparison  purposes. 

Gamma  Kerma  Factors 


21  or  36 
V”1  E;  <l>; 


i  =  1 


I  ,-i  ^Ei 


21  or  36 

I 


F.q.  1 


<l>i  AEj 


i  =  1 


Gamma  kerma  factors  for  tissue  in  a  21-group  format  were  taken  directly 
from  DLC-31  (4).  Although  DLC-31  tissue  represents  a  standard  man 
composed  of  11  elements,  the  gamma  kerma  factors  compared  reasonably 
well  with  those  calculated  from  mass  energy  transfer  coefficients  for  muscle 
(5).  Gamma  kerma  factors  for  the  elements  carbon  and  magnesium  (ion 
chamber  wall  materials)  were  also  taken  directly  from  DLC-31  (4). 
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Neutron  Kerma  Factors 


In  addition  to  gamma  kerma  factors,  selected  neutron  kerma  factors  were 
contained  in  DLC-31  (4),  but  some  materials  of  particular  interest  for  paired 
ion  chamber  constants  were  not  covered,  namely  ICRU  muscle,  muscle- 
equivalent  gas,  and  argon.  Consequently,  the  more  complete  neutron  kerma 
factor  data  of  ICRU  Report  26  (1)  were  the  basis  for  the  present  calculations. 
The  energy  grouping  of  the  ICRU-26  kerma  factor  data  did  not  match  the 
37-group  structure  of  the  neutron  spectra,  so  restructuring  of  the  groups  was 
required.  The  procedures  used  to  convert  the  ICRU-26  kerma  factor  data  to 
the  37-group  format  are  described  in  reference  6. 

Calculations 


Spectrum-averaging  calculations  for  neutron  and  gamma  kerma  factors  were 
done  with  the  AFRRI  PDP11  computer.  All  kerma  factor  data  and  all 
neutron  and  gamma  energy  spectra  were  stored  in  formatted  input  files.  A 
fortran  computer  program  named  SPT  was  written  to  complete  the  calcula¬ 
tions  and  store  the  results  in  separate  output  files.  SPT  calculated  the 
following  quantities: 

oo 

K(E)<J>(E)  dE 

0 

SAKF  =  - 

OO 

<f»(E)  dE 

0 


where  SAKF  is  the  spectrum  averaged  kerma  factor,  K(E)  and  Kj  are  the 
energy-dependent  kerma  factors,  <j>(E)  and  <pj  are  the  neutron  or  gamma 
fluence,  and  E  is  the  energy.  Note  that  the  SAKF  integrals  include  thermal 
neutrons  in  the  averaging  process.  A  listing  of  SPT  is  contained  in 
Appendix  A. 


21  or  37 

I*' 


<l>;  AEj 


i  =  1 


21  or  37 


F.q.  2 


AEj 


i  =  1 


RESULTS 


Reactor  Gamma  SAKF 


Gamma  SAKF's  for  tissue  are  shown  in  Table  2.  The  difference  between  the 
gamma  SAKF’s  for  tissue  for  the  various  reactor  configurations  can  be 
interpreted  as  changes  in  the  effective  energy  of  the  gamma  spectra.  Tissue- 
effective  energies  of  the  gamma  spectra  given  in  Table  2  represent  the  single 
photon  energies  that  would  give  the  same  tissue  kerma  factors  as  the  gamma 
spectra  in  question.  Tissue-effective  gamma  energies  ranged  from  0.6  to  1.5 
MeV  for  the  reactor  configurations,  and  the  energies  correlate  reasonably 
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well  with  the  average  energies  of  the  spectra  listed  in  Table  1.  The  fact  that 
tissue-effective  gamma  energies  were  0.1  to  0.3  MeV  lower  than  spectrum- 
averaged  energies  expresses  the  higher  attenuation  coefficients  of  lower 
energy  gamma  rays  in  tissue.  For  materials  other  than  tissue,  the  effective 
gamma  energies  could  vary  because  of  differences  in  atomic  number. 


Table  2.  Reactor  Gamma  Spectrum-Averaged  Kerma  Factors  and  Ratios* 


50 

Unshielded 

too 

Unshielded 

200 

Unshielded 

300 

Unshielded 

400 

Unshielded 

500 

Unshielded 

50 

Unshielded 

100 

Unshielded 

200 

Unshielded 

300 

Unshielded 

NA 

Pneumatic  tubes 

100 

7.5  cm  H2O 

100 

30  cm  H20 

100 

5  cm  Pb 

100 

5  cm  Pb  and  exer¬ 

cise  wheel 

100 

15  cm  Pb 

100 

15  cm  Pb  and  cave 

100 

15  cm  Pb  and  cylin¬ 

drical  phantom 

Reference  Spectra 

DLC-31 

Prompt  gamma  ray 

fission 

Cobalt-60 

,0 

4.92 

,0 

4.88 

,9 

4.47 

,8 

3.89 

,7 

3.35 

(K^-)g  units  are  10"1“  rad  •  cm2.  (K-p/Kp-)  and  (K-p'Ky]™)  are  ratios  of 
SAKF's. 


Mso  shown  in  Table  2  are  the  ratios  of  SAKF's  for  tissue  to  those  for  carbon 
and  magnesium.  The  C  and  Mg  data  are  presented  in  this  format  to  be 
directly  useful  in  dosimetry  applications.  Average  values  (±  2  standard 
deviations)  for  C  and  Mg  SAKF  ratios  for  all  reactor  configurations  are  given 
below. 

(KT/Kc)r,  =  1.1 33  ±1.0% 

F.q.  3 

(1<t/K  Mg)  r,  =  1.09  ±3.1% 
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Corresponding  kerma  factor  ratios  for  cobalt-60  gamma  rays  (1.25  MeV 
effective  energy)  are  also  shown  in  Table  2.  For  C,  cobalt-60  gamma  rays 
give  a  kerma  factor  ratio  lower  than  for  any  of  the  reactor  configurations; 
for  Mg,  it  is  higher.  This  finding  is  consistent  with  the  presence  of  Low- 
energy  photons  (below  300  keV)  in  each  of  the  gamma  spectra  of  the  reactor 
configurations.  The  average  values  (Equation  3)  for  C  and  Mg  gamma  SAKE 
ratios  differ  by  less  than  4%  from  the  corresponding  ratios  for  eobalt-60. 
However,  for  some  configurations,  the  influence  of  low-energy  gamma  rays  is 
more  pronounced  and  causes  gamma  SAKF  ratios  to  differ  from  cobalt-60 
values  by  up  to  7.6%.  This  is  the  case  at  great  distances  from  the  reactor, 
where  low-energy  gamma  rays  arising  from  the  walls  exert  an  increasing 
effect,  and  also  behind  the  15-cm  lead  shield,  where  lead  X  rays  presumably 
contribute  to  the  total  gamma-ray  spectrum.  On  the  other  hand,  the  5-cm 
lead  shield  appears  to  have  a  hardening  effect  on  the  gamma  spectrum,  as  far 
as  C  and  Mg  gamma  SAKF  ratios  are  concerned. 

Reactor  Neutron  SAKF  for  Tissue 

Neutron  spectra  calculated  from  3-D  models  are  more  reliable  than  those 
obtained  from  1-D  models.  Consequently,  the  present  calculations  used  3-D 
neutron  spectra  for  all  reactor  configurations  for  which  they  were  available 
(see  Table  1).  Comparison  of  SAKF's  for  ICRTJ  muscle  obtained  from  1-D  and 
3-D  spectra  is  shown  in  Table  3  for  those  configurations  for  which  both  were 
available.  For  noncomplicated  configurations  such  as  unshielded  at  a 
distance  of  100  cm,  close  agreement  occurs  between  the  1-D  and  3-D  results. 
This  agreement  confirms  the  adequacy  of  the  1-D  spectra  for  SAKF  calcula¬ 
tions  at  the  other  unshielded  distances.  Likewise,  the  reasonable  agreement 
between  the  1-D  and  3-D  results  for  the  reactor  shielded  bv  30  cm  of  water 
justifies  the  use  of  a  1-D  spectrum  for  the  7. 5-cm  water  shield.  On  the  other 
hand,  the  lack  of  a  3-D  spectrum  for  the  pneumatic  tubes  leaves  open  the 
question  of  the  adequacy  of  the  1-D  neutron  spectrum  at  that  location. 

Table  3  also  presents  SAKF  values  for  ICRU  muscle  for  the  five  measured 
neutron  spectra.  These  values  differ  by  2%-58%  from  those  determined  from 
the  3-D  calculated  spectra.  The  closest  agreement  is  for  the  lead  cave 
configuration,  for  which  the  calculated  and  measured  neutron  spectra  plotted 
in  reference  3  are  nearly  coincident.  For  the  worst  case,  the  cylindrical 
phantom,  the  spectrum  plots  show  a  clear  excess  of  neutrons  below  300  keV 
in  the  calculated  spectrum.  Reference  3  explains  this  as  being  due  to 
shortcomings  in  the  calculated  spectrum,  on  the  assumption  that  the 
measured  neutron  spectrum  is  the  more  accurate.  Thus,  the  SAKF  differ¬ 
ences  between  the  3-D  and  measured  spectra  can  be  interpreted  as  an 
estimate  of  the  uncertainty  in  SAKF  values  due  to  uncertainties  in  the 
neutron  spectra. 


Table  3.  1CRU  Muscle  SAKF  for  Calculated  and  Measured 
Reactor  Neutron  Spectra 


Table  4  presents  the  reactor  neutron  SAKF's  for  three  tissuelike  materials. 
ICRU  muscle,  the  generally  accepted  reference  material  for  neutron  dosime¬ 
try  (7,8)  is  used  as  such  here.  ICRU  tissue  differs  from  ICRU  muscle  by  a  2% 
lower  hydrogen  content  (by  weight),  and  the  calculated  kerma  factors  reflect 
this  difference.  DLC-31  tissue  represents  standard  man  composed  of  11  ele¬ 
ments,  but  the  neutron  SAKF  values  differ  negligibly  from  those  for  ICRU 
muscle  or  ICRU  tissue.  A-150  is  the  tissue-equivalent  plastic  used  for  tissue- 
equivalent  ion  chamber  construction. 

Tissue-effective  energies  of  the  neutron  spectra  tabulated  in  Table  4  were 
determined  in  the  same  way  as  those  of  the  gamma  spectra.  That  is.  tissue- 
effective  spectrum  energies  were  taken  as  the  monoenergetic  neutron 
energies  giving  the  same  tissue  kerma  factors  as  the  spectra  in  question.  The 
tissue-effective  neutron  energies  display  a  decreasing  trend  with  increasing 
free-in-air  distance  from  the  reactor  (due  to  the  increased  importance  of 
scattered  neutrons),  and  an  increasing  trend  for  water-shielded  configurations 
(due  to  a  beam-hardening  effect).  Two  phenomena  are  noteworthy:  (a)  The 
tissue-effective  neutron  energies  are  not  equal  to.  nor  do  they  correlate  very 
well  with,  the  average  neutron  energies  given  in  Table  1.  This  finding 
emphasizes  the  fact  that  neutron  spectra  of  complex  shape  cannot  be 
represented  well  by  any  single  energy  descriptor,  (b)  The  tissue-effective 
neutron  energies  are  low  at  the  two  points  with  hydrogenous  surrounds,  i.e.. 
pneumatic  tubes  and  cylindrical  phantom.  Ponsideration  of  the  shapes  of  the 
two  neutron  spectra  plotted  in  Figure  1  clarifies  this  phenomenon.  The 
spectra  pertain  to  the  reactor  shielded  with  15  cm  of  Pb,  at  a  distance  of  100 
cm.  One  spectrum  represents  the  case  with  no  phantom  present:  the  other 
spectrum  is  that  at  the  center  of  the  1 8-cm-diameter  cylindrical  phantom. 
Within  the  phantom  the  total  neutron  flux  below  0.01  WeV  is  8.8  times  the 
flux  without  the  phantom,  while  the  flux  above  0.01  YleV  is  decreased  by  a 
factor  of  5.  It  is  the  large  low-cncrgv  component  of  the  neutron  spectrum 
within  the  phantom  that  causes  such  low  tissue-effective  neutron  energy. 


Table  4.  Reactor  Neutron  Spectrum-Averaged  Kerma  Factors  for  Tissue 


(from  reference  1) 


Reactor  Neutron  SAKF  Ratios 


Distance 
to  Core 
Room  (cm) 

Configuration 

SAKF  (10' 

10  rad  cm2) 

Cl 

/J 

MeVeff 

Muscle^ 

(ICRU) 

Tissue  4. 
(ICRU)  1 

Tissue 

(DLC-31)t 

A-150  . 
Plastic  3 

1  50 

Unshielded 

0.8 

19.8 

19.6 

19.6 

20.0 

too 

Unshielded 

0.7 

20.3 

20.2 

20.0 

20.6 

200 

Unshielded 

0.6 

17.6 

17.4 

17.4 

17.8 

300 

Unshielded 

0.4 

15.4 

15.2 

15.2 

15.6 

'  * 

400 

Unshielded 

0.4 

13.6 

13.4 

13.4 

13.7 

500 

Unshielded 

0.3 

12.1 

11.9 

11.9 

12.2 

2  50 

Unshielded 

0.4 

14.2 

14.0 

14.0 

14.4 

too 

Unshielded 

0.4 

14.4 

14.2 

14.2 

14.5 

.  • 

200 

Unshielded 

0.3 

12.1 

11.9 

1 1.9 

12.2 

300 

Unshielded 

0.2 

10.6 

10.5 

10.5 

10.7 

s 

NA  NA 

Pneumatic  tubes 

0.04 

3.83 

3.70 

3.75 

3.37 

•  « 

m 

1  too 

7.5  cm  H2O 

1.0 

23.7 

23.4 

23.4 

24.1 

-*\ 

too 

30  cm  H2O 

2.7 

34.4 

33.9 

34.0 

35.1 

too 

5  cm  Pb 

0.8 

20.5 

20.2 

20.2 

20.4 

J 

too 

5  cm  Pb  and  exer- 

r- 

cise  wheel 

0.3 

12.0 

11.8 

11.8 

12.0 

1 

too 

15  cm  Pb 

0.6 

17.0 

16.8 

16.8 

16.8 

too 

15  cm  Pb  and  cave 

0.4 

13.6 

13.4 

13.4 

13.5 

too 

15  cm  Pb  and  cylin- 

•  • 

drical  phantom 

0.03 

2.85 

2.81 

2.77 

2.83 

•J 

Reference  Spectra 

DLC-31  fission 

0.8 

20.4 

20.1 

20.1 

20.4 

tfl 

DLC-31  thermonuclear 

0.6 

18.1 

17.8 

17.8 

18.3 

’  « 

ENDF  fission 

1.16 

28.3 

27.9 

27.9 

28.4 

DLC-31  14  MeV 

14. 

64.6 

63.9 

63.2 

67.4 

0m 

*  Muscle  (ICRU):  10.2 %  H,  12.3%  C,  3.5%  N. 

72.9%  O. 

1.1%  (Na  * 

Mg  - 

m 

P  +  S  +  K  +  Ca)  (from  reference  1) 

1 

4- 

4 

Tissue  (ICRU): 

10.0%  H,  14.9%  C,  3.5%  N, 

71.6%  O  (from  reference  l) 

Tissue  (DLC-31):  10%  H,  24%  C, 

2.9%  N,  60%  O,  0.20%  Na.  0.30%  Mg. 

1.1%  P,  0.24%  S 

0.20%  K, 

1.2%  Ca, 

0.20%  Cl 

(from 

X 

reference  5  of  reference  4) 

■ 

§ 

A-150  plastic: 

10.1%  H,  77.6% 

C,  3.5%  N 

5.2%  O. 

1.8%  Ca.  1 

.7%  F 

«* 

SAKE  results  for  the  three  ion  chamber  gas  materials  are  shown  in  Table  5. 
For  TE  gas,  all  the  reactor  configurations  can  be  represented  by  a  single 
(Kj/K'pp;  gas^N  value  of  0.983  ±  0.5%  (2  SD).  For  CO2  and  Ar,  the  SAKF 
ratios  can  be  summarized  as  9.78  ±  10%  and  70.7  ±  27%.  respectively.  The 
variability  among  the  CO2  and  Ar  SAKF  ratios  indicates  that  individually 
tabulated  values  should  be  used  for  each  reactor  configuration. 
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Table  5.  Reactor  Neutron  SAKF  Ratios 


Room 

Distance 
to  Core 
(cm) 

Configuration 

(Kt/Kte  sas)N* 

(KT/Kco2> 

(Ky/  K 

1 

50 

Unshielded 

0.980 

9.57 

66.0 

100 

Unshielded 

0.979 

9.50 

64.3 

200 

Unshielded 

0.983 

9.62 

66.7 

300 

Unshielded 

0.981 

9.68 

67.5 

400 

Unshielded 

0.986 

9.78 

69.0 

500 

Unshielded 

0.984 

9.84 

70.3 

2 

50 

Unshielded 

0.979 

9.59 

66.0 

100 

Unshielded 

0.986 

9.66 

65.8 

200 

Unshielded 

0.984 

9.68 

67.2 

300 

Unshielded 

0.981 

9.81 

69.3 

NA 

NA 

Pneumatic  tubes 

0.982 

10.21 

73.2 

1 

100 

7.5  cm  HoO 

0.979 

9.26 

61.2 

100 

30  cm  H<>0 

0.983 

8.49 

5  5 . 5 

100 

5  cm  Pb 

0.986 

9.67 

69.7 

100 

5  cm  Pb  and  exer- 

cise  wheel 

0.984 

9.92 

73.6 

100 

15  cm  Pb 

0.983 

10.18 

83.7 

100 

15  cm  Pb  and  cave 

0.986 

10.38 

92.5 

100 

15  cm  Pb  and  evlin 

- 

drical  phantom 

0.986 

11.00 

91.1 

Reference  Spectra 

j 

DLC-31  fission 

0.986 

9.86 

76.7 

DLC-31  thermonuclear 

0.989 

6.24 

50.7 

ENDF  fission 

0.983 

9.53 

65.2 

DLC-31  14  MeV 

0.977 

3.40 

26.5 

*  TE  gas:  10.2%  H,  45.6%  C,  3.5%  N.  40.7%  O 


The  constancy  of  (Kt/Kte  gas^N  among  the  reactor  configurations  was  to  be 
expected,  because  this  ratio  varied  only  from  0.95  to  1.01  over  the  entire 
range  of  37-energv  groups.  For  the  other  neutron  SAKF  ratios,  uncertainties 
can  be  estimated  by  again  comparing  results  from  calculated  and  measured 
spectra,  as  shown  in  Table  6.  Although  individual  SAKF  values  have 
uncertainties  of  2%-58%,  SAKF  ratios  determined  from  calculated  and 
measured  spectra  differ  by  only  l%-8%. 


Table  6.  SAKF  Ratios  for  Calculated  and  Measured  Neutron  Spectra 


Room 

Distance 
to  Core 
(cm) 

Configuration 

Neutron  Spectrum 

3-D  MEAS 

Differ¬ 
ence  (%) 

Kj/K  Ar 

1 

100 

Unshielded 

64.3 

67.5 

-5 

2 

100 

Unshielded 

65.8 

67.7 

-3 

1 

100 

5  cm  Pb  and  exer- 

cise  wheel 

73.6 

74.6 

-1 

l 

100 

15  cm  Pb  and  cave 

92.5 

91.4 

-1 

1 

100 

15  cm  Pb  and  cylin- 

drical  phantom 

91.1 

84.3 

*8 

kT/kCO? 

1 

100 

Unshielded 

9.50 

9.66 

_o 

2 

100 

Unshielded 

9.66 

9.73 

-1 

1 

100 

5  cm  Pb  and  exer- 

cise  wheel 

9.92 

10.00 

-1 

1 

100 

15  cm  Pb  and  cave 

10.38 

10.30 

+1 

1 

100 

15  cm  Pb  and  cylin- 

drical  phantom 

11.0 

10.56 

*4 

DISCUSSION 


The  gamma  SAKF  ratios  for  C  and  Mg  (ion  chamber  wall  materials)  vary  only 
slightly  among  the  reactor  configurations,  and  average  values  represent  these 
quantities  with  fairly  high  precision.  Historically  it  has  been  questioned 
(9,10)  whether  the  AFRRI  reactor  gamma  spectra  can  be  approximated  by 
cobalt-60  radiation  (1.25  MeV).  The  present  calculations  confirm  that  this  is 
nearly  the  case,  since  the  average  C  and  Mg  SAKF  ratios  differ  less  than  4% 
from  the  values  applicable  to  cobalt-60.  However,  in  some  configurations,  a 
low-energy  photon  component  in  the  reactor  gamma  spectra  does  cause  C  and 
Mg  SAKF  ratios  to  differ  by  up  to  7.6%  from  corresponding  values  for 
cobalt-60. 

The  constancy  of  the  neutron  SAKF  ratio  for  TE  gas  between  the  various 
reactor  configurations  was  to  be  expected,  because  the  compositions  of  TE 
gas  and  ICRU  muscle  match  closely.  For  CO2  and  Ar,  the  neutron  SAKF 
ratios  are  nearly  constant  for  all  free-in-air  unshielded  reactor  configura¬ 
tions,  but  the  ratios  differ  considerably  (5%-40%)  for  more  complex  reactor 


configurations. 


The  d  fferences  between  SAKF  values  derived  from  3-D  calculated  neutron 
spectra  and  those  derived  from  measured  neutron  spectra  allow  estimation  of 
the  uncertainties  in  the  reported  SAKF  values  arising  from  spectrum  uncer¬ 
tainties.  For  ICRU  muscle,  the  observed  2%-58%  differences  place  a  limit  on 
dosimetric  accuracy  that  can  be  obtained  using  a  direct  fluenee-to-dose 
method.  This  emphasizes  the  importance  of  using  ionization  chambers  for 
dose  determination  since  the  overall  uncertainty  of  those  chambers  is  about 
10%  (8).  In  this  regard,  the  SAKF  ratios  for  CO2  and  Ar  differ  by  only  1%- 
8%  between  the  3-D  calculated  and  measured  neutron  spectra.  The  improved 
precision  for  SAKF  ratios  over  individual  SAKF  values  is  due  to  the  similarity 
in  energy  response  of  the  kerma  factors  for  the  different  materials.  This 
improved  precision  is  indeed  fortunate,  because  the  SAKF  ratios  are  required 
to  evaluate  paired  ionization  chamber  response  constants  for  reactor  neutron 
and  gamma  dosimetry. 

Considering  the  SAKF  results  for  the  reference  spectra,  it  was  apparent  that 
for  ICRU  muscle  and  TE  gas,  the  prompt  gamma-ray  fission  source  as  well  as 
the  fission  and  thermonuclear  sources  fell  within  the  range  of  values  obtained 
for  AFRRI  reactor  configurations.  In  particular,  the  15-cm  lead  shield  gave 
the  same  tissue-effective  neutron  energy  as  the  thermonuclear  source,  and 
the  5-cm  lead  shield  matched  the  fission  source.  However,  these  apparent 
similarities  should  be  viewed  with  caution  because,  as  pointed  out  above,  no 
single  energy  descriptor  can  adequately  represent  a  complex  energy 
spectrum. 

The  results  of  this  study  demonstrate  the  need  for  further  neutron  spectrum 
calculations  or  measurements  for  AFRRI  reactor  configurations  at  depth  in 
tissue-equivalent  or  anatomical  phantoms.  The  large  change  in  the  neutron 
spectrum  within  the  18-cm-diameter  phantom  shows  the  importance  of  also 
studying  phantoms  of  different  sizes  and  shapes.  The  very  low  tissue- 
effective  energy  of  the  neutron  spectrum  within  the  phantom  also  suggests 
the  need  for  more  sophisticated  dosimetric  tools  (such  as  tissue-equivalent 
proportional  counters)  to  supplement  the  paired  ionization  chambers  and  to 
better  quantitate  the  dose  deposition  process. 
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LATEST  CORRECTION  17  AUG  82 

THIS  PROCRAH  WILL  COMPUTE  TWO  SPECTRUM  DEPENEENT  FUNCTIONS: 

A  SPECTRUM  WEICHTEB  FUNCTION  AND  A  SPECTRUM  FOLDED  FUNCTION, 

AND  THE  ESTIMATED  VARIATION  OF  EACH 

AUTHOR:  KEI  FERLIC  CIRCA  1981  OR  82 

THIS  PROGRAM  IS  FOR  EVALUATING  THE  INTEGRALS  (BY  SUMMATION) 
NECESSARY  FOR  SPECTRUM  AVERAGED  QUANTITIES 

NOTES 

1.  SPECTRUM  INPUT  DATA  CAN  BE  A  FROM  AND  BACK  SPECTRUM 

2.  SPECTRUM  INPUT  DATA  CAN  BE  LETHARGY,  PHI(E)DE  OR  PHI(E) 

3.  THE  DELTA  SPECTRUM  IS  THE  FACTIONAL  CHANGE  EXPECTED  IN  NEUTRON 
POPULATION  FOR  A  PARTICULAR  ENERCY  GROUP 

4.  THE  DELTA  FUNCTION  IS  THE  FRACTIONAL  CHANGE  EXPECTED  IN  THE 
FUNCTION  PER  ENERGY  GROUP 

5.  ENERGY  GROUP  BOUNDS  ARE  THE  UPPER  EOUNDS 

6.  ENERGY  GROUP  1  IS  THE  MAXINUM  ENERCY  GROUP.  I.E.  GROUP  1  OFFER 
BOUND  IS  19-6  MEV 


HIGHEST  CONTROL  STATEMENT  IS  238 
HICHEST  FORMAT  STATEMENT  IS  1217 


BYTE  FILE1 (3C)  1FILE1 
BYTE  F I LE2 ( 30 )  'FILE2 
BYTE  FILE3(3C)  'FILE3 
BYTE  F I LE4 ( 30 )  1FILE4 
BYTE  FILE5(3C)  1FILE5 
BYTE  F I LE6 ( 30 )  1FILE6 
DIMENSION  PHITC38) 
DIMENSION  PHIF(38) 
DIMENSION  PHIBC38) 
DIMENSION  ICR0UP(38) 
INTEGER  NROW 
INTEGER  ICOL 
INTECER  ICNT 
INTEGER  NSPT 
INTEGER  NTYPE 
DIMENSION  ENCP(38) 
DIMENSION  FNCT ( 38 ) 
DIMENSION  DE(38) 
DIMENSION  DU ( 38 ) 

REAL  DUBWI 
REAL  DLBWI 
REAL  DUBNI 
REAL  DLBNI 
INTEGER  IUBGWI 
REAL  UBFWI 
INTEGER  IUBFGW 
INTEGER  IUBGNI 
REAL  UBFNI 
INTEGER  IUBFGN 
INTEGER  ILECWI 
REAL  LBFWI 
INTEGER  ILEFGW 
INTEGER  ILBFGN 
REAL  DEFNI 


■FLUX  SPECTRUM  FILE 
■ENERGY  GROUP  FILE,  UFPER  BOUNDS 
■FUNCTION  FILE  EY  CROUP 
■OUTPUT  FILE 
■DELTA  SPECTRUM  FILE 
■DELTA  FUNCTION  FILE 
! TOTAL  INPUT  FLUX 
!  TOTAL  FRONT  FLUX 
•TOTAL  EACK  FLUX 
•CROUP  NUMBER 

I  NUMBER  CF  GROUPS  NUMBER  CF  INPUT  ROWS 
INUMBER  CF  COLUMNS  IN  INPUT  DATA 
! COUNTER  CF  FCR  FILE  NAMES 
1TYPE  OF  FLUX  INPUT  IE  ENERCY, LETHARCY 
1TYPE  OF  INPUT  SPECTRUM  TCTAL  VS  FRONT  ♦  F 
! ENERGY  GROUP  UPPER  BOUNDS 
•FUNCTION  GROUPS 
! CHANGE  IN  ENERCY 
! CHANCE  IN  LETHARGY 

1  DESIRED  UPFER  BOUND  WEIGHTING  INTEGRAL 
1  DESIRES  LOWER  BOUND  VEICHTINC  INTEGRAL 
1  DESIRED  UPPER  BOUND  NORMALISING  INTFCRAL 
! DESIRED  LOWER  FOUND  NORMALIZING  INTFCRAL 
•UPPER  POUND  GROUP  WEIGHTING  INTEGRAL 
•UPPER  BOUND  FRACTION  WEIGHTING  INTEGRAL 
•UPPER  BOUND  FRACTIONAL  GROUP  WEICHTINC  IN 
•UPPER  POUND  CROUP  NORMALIZING  INTECRAL 
! UPPER  BOUND  FRACTION  NORMALIZING  INTECRAL 
1  UPPER  BOUND  FRACTIONAL  CF.CUF  NORMALIZING 
•LOWER  BOUND  GROUP  WEICHTINC  INTECRAL 
•LOWER  FOUND  FRACTION  WEICHTINC  INTECRAL 
•LOWER  BOUND  FRACTIONAL  CROUP  WEICHTINC  IN 
1  LOWER  BOUND  FRACTIONAL  GROUP  NORMALIZING 
•LOWER  POUND  FRACTION  NORMALIZING  INTEGRAL 


v. 


INTEGER  ILBGNI 
SEAL  HIST 
HEAL  SINT 
SEAL  FNCEFF 
ISTEGEE  NS0H1 
INTECEH  R0V2 
INTEGER  IH0V3 
INTEGER  NH0H5 
INTECES  NH0H6 
DIMENSION  FRCTW { 58 ) 
DIMENSION  FRCTR ( 38 ) 
DIMENSION  DELTA (38) 
REAL  UTAIL 
BEAL  LTAIL 
INTEGER  SAF 
INTEGER  VSAF 
INTEGER  SFF 
INTEGER  VSFF 
INTEGER  IKPUT1 
INTECER  INPUT2 
INTECER  IN PUT 3 
DIMENSION  DFF ( 38 ) 
DIMENSION  DSF(38) 
DIMENSION  PHIENR(38) 
REAL  UBFEH 
REAL  UBFEN 
REAL  LBFEW 
REAL  LBFEN 
DIMENSION  DEW (38) 
DIMENSION  DEN (38) 
DIMENSION  PNPHI ( 38) 
DIMENSION  PWFNT(3S) 
DIMENSION  PHSUM(38) 
REAL  VFFNCT 
DIMENSION  TERM' (38) 
DIMENSION  TERM2(38) 
DIMENSION  TERM3(38) 
DIMENSION  PWFPHI(38) 
DIMENSION  PWEFNT ( 38 ) 
DEAL  VVFNT 


'LOWER  BOUND  CROUP  NORMALIZING  INTEGRAL 
'INTECRAL  OF  THE  WEIGHTING  INTERCRAL 
'INTEGRAL  OF  THE  NORMALIZING  INTEGRAL 
'EFFECTIVE  VALUE  OF  FUNCTION 
'ROWS  IN  INPUT  SPECTRUM  FILE 
'ROWS  IN  INPUT  ENEGRY  GROUP  FILE 
'ROWS  IN  INPUT  FUNCTION  CROUP  FILE 
'ROWS  IN  INPUT  DELTA  SPECTRUM  FILE 
'ROWS  IN  INPUT  DELTA  FUNCTION  FILE 
'FRACTION  OF  WEIGHTING  INTEGRAL  FOR  GROUP  I 
'FRACTION  OF  NORMALIZING  INTECRAL  FOR  GROUP  I 
!DE  OR  DU  IN  INTEGRAL  FOR  GROUP  I 
'UPPER  TAIL  OF  INTECRAL  BETWEEN  GROUPS 
'LOWER  TAIL  OF  INTEGRAL  BETWEEN  CROUPS 
'CONTROL  FOR  SPECTRUM  AVERAGED  FUNCTIONS 
'CONTROL  FOR  VARIATION  OF  SPECTRUM  AVERACED 
'CONTROL  FOR  SFECTRUM  FOLDED  FUNCTION 
'CONTROL  FOR  VARIATION  OF  SPECTRUM  FOLDED 
'DELTA  (VARIATION)  FILES  CONTROL  CHARACTER 
'BOUNDS  ON  NORMALIZING  INTEGRAL  CONTROL 
'TEMPORARY  NORMALIZING  INTEGRAL  CONTROL 
'DELTA  FUNCTION  ENTRY 
'DELTA  SPECTRUM  ENTRY 
! PH I ( E )  FOR  EACH  GROUP 
'FRACTION  OF  ENERCY  UPPER  END  WCT  I  NT 
'FRACTION  OF  ENERGY  UPPER  END  NCR  INT 
'FRACTION  OF  ENERGY  LOWER  END  WCT  INT 
'FRACTION  OF  ENERGY  LOWER  END  NOR  INT 
'DELTA  E  FOR  WCT  INT  IN  VARIATION  CALCULATION 
'DELAT  E  FOR  NCR  INT  IN  VARIATION  CALCULATION 
'PARTIAL  WCT  WITH  RESPECT  TO  PHI 
'PARTIAL  WCT  WITH  RESPECT  TO  THE  FUNCTION 
'SUM  OF  PART I A LS  WITH  RFSP  TO  PHI  AND  FUNCT 
’VARIATION  OF  FOLCEl  FUNCTION 
’FIRST  TERM  IN  CALCULATION  OF  VAR  OF  WCT  FNCT 
’SECOND  TERM  IN  CALCUL  OF  VAR  OF  WGT  FNCT 
'THIRD  TERM  IN  CALCUL  OF  VAR  OF  WCT  FNCT 
'PARTIAL  OF  WCT  FUNCT  WITH  RESPECT  TO  PHI  SOFT 
'PARTIAL  OF  WGT  FNCT  WITH  RESPECT  TO  FNT  SCRI 
'VARIATION  OF  WEICHTED  FUNCTION 


C 

C  COMMON  STATEMENTS  FOR  SUBROUTINES  UPBND  AND  LWBND 

C 

COMMON  IXCP,XDUB,XE(38) , IXUBG , IXUBFG , XUBF 
COMMON  XDLB, IXLBG , IXLBFG , XLBF , XUBFE , XLEFE 

C 

0  DETERMINATION  OF  CALCULATIONS  TO  BE  PERFORMED 

C 

TYPE  120C 

1200  FORMATC  DO  YOU  WANT  A  SPECTRUM  AVERACED  FUNCTION”  YES-1  .NO-C. 

1  OR  OR’) 

ACCEPT  1201, SAF 

1201  FORMAT ( 13 ) 

IF (SAF  .EG.  C)CC  TO  200 
TYPE  1202 

1202  FORMATC  DC  YOU  WANT  THE  ESTIMATE!  VARIATION”  YES-1,  NC-C  CP  CR 1  ' 
ACCEFT  1201, VSAF 

20C  CONTINUE 


16 


TYPE  1203 

1203  FORMATC  DO  YOU  WANT  A  SPECTRUM  FOLDED  FUNCTION’  YES-1  ,N0»0  OR  OR') 
ACCEPT  1201 ,SFF 

IF(SFP  .EC.  0)C0  TO  201 
TYPE  1204 

1204  FORMATC  DO  YOU  WANT  THE  ESTIMATED  VARIATION’  YES-1,  NO-C  OR  CR  * ) 
ACCEPT  1201.VSFF 

2C1  CONTINUE 

IF (SAP  .NE.  SFF }C0  TO  203 
IF(SAF  .EC.  0)C0  TO  2C2 
GO  TO  203 
202  TYPE  1205 

1205  FORMATC  WHAT  DO  YOU  V ANT’?’  BYE,  BYE1) 

CO  TO  140 

2C3  CONTINUE 

C 

C  PROGRAM  CONTROL  FOR  CALCULATIONS,  DETERMINATION  OF  ACTUAL  INPUT  DATA 

C 

C 

C  INPUT1  DELTA  FILES  INPUT, 

C  INPUT2  BOUNDS  ON  NORMALIZING  INTECRAL  INPUT 

C 

IF(SAF  .EC.  C )C0  TO  2C4 
INPUT  1-0 
INPUT2- 1 

IFfVSAF  .EC.C)CO  TO  2C4 
INPUT1 -1 
2C4  CONTINUE 

IF(SFF  .EC.  G.'CC  TO  205 

INPUT1-C 

INPUT3-0 

IF(VSFF  .EC.O.'GO  TO  2C5 
1NPUT1-1 
2C5  CONTINUE 

IF( INPUT2  .EC.  INPUT3)CC  TO  229 
IKPUT2-1 
229  CONTINUE 

C 

C  IDENTIFY  NUMEER  OF  GROUPS 

C 

TYPE  1116 

1116  FORMATC  WHAT  IS  THE  NUMBER  CF  GROUPS  IN  SPECTRUM: <13 > ‘ '/ 

ACCEPT  1006, NROW 

C 

C  SET  LOW  BOUND  ON  LOWEST  ENERGY  GROUP 

C 


Ir (NROW  .EC.  37 )C0  TO  119 
IF (NROW  .EC.  21 }G0  TO  12C 
TYPE  1120 

1120  FORMATC  WHAT  IS  LOWER  BOUND  (IN  EV )  CN  THE  LOWEST  ENEECY' 
1 '  GR0UP?<E10.4>‘  ) 

ACCEPT  1121,  ENGP(NROW*C 

1121  FORMAT(EIO.e) 

GO  TO  121 

119  ENGP( 38)-1 .CE-5  'IN  EV 

CO  TO  121 


4 

'r 
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120  £NCP(22)*1 .0E-4  'IN  EV 

121  CONTINUE 

C 

C  FILE!  NAME:  SPECTRUM  FILE 

C 

TYPE  1000 

1000  FORMATC  ENERCY  SPECTRUM  FILE  NAME  <FILENAMEX.TYP> :  ') 

ACCEPT  1001 ,ICNT, FILE! 

10C1  FORM AT ( 0 , 30A 1 ) 

FILE1 ( ICNT*1 )*0 
TYPE  1005 

1005  FORMATC  SPECTRUM  INPUT:  FRONT  ♦  BACK  -  C,  TOTAL  -  1  ' ) 
ACCEPT  1C06,  NTYPE 

1006  F0RMAT(I3 ) 

TYPE  1115 

1115  FORMATC  IS  SEPCTRUM  INPUT:  PHI (l')»1 ,  PHI  ( E)DE-C,PHI I £)  —  1  ' ) 
ACCEPT  10C6.NSPT 


FILE2  NAME:  ENERCY  CROUP  UPPER  EOUNE  FILE 


TYPE  1CC3 

1CC3  FORMATC  ENTER  ENERCY  GROUP  (UFPER  ECUNP}  FILE:  <FILENAMEX.TYP> 1 ) 
ACCEPT  1CC1 .ICNT.FILE2 
FILE2(ICNT*1 )*C 

C 

C  FILE3  NAME:  FUNCTION  CROUP  FILE 

C 

TYPE  1GC4 

1CC4  FORMATC  ENTER  FUNCTION  CROUP  FILE:  <’FILENAMEX.TYP>  ' } 

ACCEPT  1 CC 1 , ICNT.FILE3 
FILE3( ICNT’1 )*C 

C 

C  FILE4  NAME:  OUTPUT  FILE  NAME 

C 

TYPE  1196 

1196  FORMATC  ENTER  OUTPUT  FILE  NAME:  <FILENAMEX.TYP>  C 
ACCEPT  1001 ,INCT,FILE4 
FILE4 ( INCT*1 )-0 

C 

C  FILE5  NAME:  DELTA  SPECTRUM  FILE 

C 

IF( INPUT1  .EC.  C)CO  TO  2C6 
TYPE  1 2C6 

1206  FORMATC  ENTER  DELTA  SPECTRUM  FILE  NAME  •  FI LENAMEX . TYP> :  ' ) 

ACCEPT  10C1 .INCT.FILE5 
FILE5(INCT*1 )-C 

C 

C  FILE6  NAME:  DELTA  FUNCTION  FILE 

C 

TYFE  1207 
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1207  formatC  enter  delta  function  file  name  <filenamex.typ> 
ACCEPT  1CC1 ,INCT,FILE6 
FIL£6(INCT*1 )*0 
2C6  CONTINUE 

C 

C  EXTRACT  DATA  FROM  SPECTRUM  FILE;  FILE1 

C 

0PEN(UNIT*1 .NAME-FILE1 .TYPE* ‘OLD ’ ) 

READd  ,1007)  IC0L.NR0V1 

1007  F0R«AT(2I4) 

IF(NR0N1  .NE.  NR0W)G0  TO  112 
IF(NTYPE  .EC.  C)GO  TO  100 
DO  101  I-1.NR0W 
READd, 1008)  IGROUP( I).PHIT(I) 

1008  F0RMAT(I3,E1 2.4 ) ' 

101  CONTINUE 
GO  TO  1C2 

ICO  DO  103  I-1.NRCW 

READO  ,  1 CC9 )ICRCl'F ( I ) , PKIF( I ) ,  PKIB(  I ) 

1009  FCRMAT( I3.2E1 2. 4 ) 

1 03  CONTINUE 

102  CONTINUE 
CLCSE(UNIT*1  ) 

C 

C  EXTRACT  DATA  FROM  ENERGY  GROUP  FILE 

C 

0PEN(UNIT-1 .NAME-FILE2.TYPE* 'OLD' ' 

READ(  1  , 1007 )  ICO  I..NRCW2 
IF(NR0W2  .NE.  NRCV'CO  TO  114 
DO  104  I-1.NR0V 

READ(1,1CC6)  ICRCUPCD.ENCPa) 

1C4  CONTINUE 

C LOSE (UN  IT* 1 ) 

C 

C  EXTRACT  LATA  FROM  FUNCTION  FILE: FILE? 

C 

OPEN ( UNIT* 1 ,NAME-FILE3,TYPE*'0LD’ ) 

READ( 1 , 1007), IC0L.NRCN3 
IF ( NR0W3  .NE.  NROV )C0  TO  119 
DO  105  1*1 ,  NROV 

READ( 1 ,1008),  ICROUP( I ) , FNCT ( I ) 

105  CONTINUE 

CL0SE(UNIT*1 ) 

C 

C  EXTRACT  LATA  FROM  DELTA  SPECTRUM  FILE:  FILE1 

C 

IF f INFUT1  .EC.  C )GC  TO  207 
0PEN(UNIT-1 .NAME-EILE5. TYPE- ’ILL ' ■ 

READ(  1 , 1 007 ) ,  ICCL,  NF.t  V5 
IF(NR0W5  .NE.  NPCh’GC  TO  DO¬ 
DO  2C9  1-1  ,NRO> 

READ(1  ,  1CCe),IGRCUPd),DSFd 
209  CONTINUE 

ULCSEf  UNIT-1  : 
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c 

C  EXTRACT  DATA  FROM  DELTA  FUNCTION  FILE:FILE6 

C 

OPEN(UNIT-1 .NAME-FILE6, TYPE- ‘OLE  C' 

READ(1 , 1007) , IC0L.NR0V6 
IF(NR0W6  .NE.  »ROM)GO  TO  21C 
DO  211  1-1 ,NROV 
READ ( 1 , 1008) , ICROUP(I ),EFF( I ) 

211  CONTINUE 

CLOSE(UNIT- 1 ) 

2C7  CONTINUE 

GO  TC  141 

C 

C  CHECK  FOR  PROPER  MUMPER  OF  GROUPS  IN  EACH  INFUT  FILF;  NEGATIVE  ANSWFF 

C 

112  TYPE  1 1 17.NR0V1 
GC  TO  113 

114  TYPE  1117.NR0W2 
GC  TO  113 

115  TYPE  1 1 1 7.NR0W3 
GC  TO  113 

2l&  TYFE  1117.NRCW5 

GC  TC  113 

21C  TYPE  1 1 1 7, NR0W6 

1117  FORMATC  ACTUAL  GROUPS  IN  FILE ‘ , 13  ) 

113  TYPE  1 1 99.NR0N 

1 1  V9  FORMATC  DOES  NCT  EC  UAL  IDENTIFIED  SFECTR'.'M  GRCUFS:  ' ,  13! 

CLCSE(UNIT-1 ) 

GO  TO  14C 

C 

C  CALCULATE  TOTAL  FLUX  IF  FRONT  AND  BACK 


141  IF(NTYFE  .EC.  1 )G0  TO  106 

DO  1C7  1-1,  NRGW 
PHIT( I )-FHIF(I )*PHIB( I) 
107  CONTINUE 

106  CONTINUE 


C 

C  CALCULATE  DU  OR  DE  FOR  THE  SPECTRUM 

C 

C  CLD  LINE  IF(NSPT)  116,117.116 

C 

C  CHANGE  IN  LETHARGY/ENERGY 

C 

116  DC  1C&  1*1 , NFOW 

DU  C I ; -LCC ( ENGFf I ' /EHGF( 1*1 )) 

106  CONTINUE 

C  OLE  LINE  GC  TO  li¬ 
nt  EC  1C9  1-1 , NRCW 

EE(I)*ENGF(I )-ENCP( I ♦'  , 

ICS  CONTINUE 

117  CONTINUE 


•  C 

cl 


c 

c 

c 


ACCEPT  BOUNDS  ON  WEICHTINC  AND  NORMALIZING  INTEGRAL 


TYPE  1010 

1010  FORKATC  WHAT  IS  UPPER  BOUND  OF  WEIGHTING  INTECRAL  IN  EV:' 

1  ' <E10.4> 1 ) 

TYPE  1122 

1122  FORMATC  NEUTRON  37  GROUP  MAX  IS  19.6E-6  EV;  PHOTON  21  CROUP  MAX' 
1  'IS  14.0E»6  EV) 

ACCEPT  1011.DUBWI 
IC11  FORMAT! El  0. 4 ) 

TYPE  1012 

1012  FORMATC  WHAT  IS  LOWER  BOUND  OF  WEIGHTING  INTEGRAL  IN  EV:  1 
1  ' <E1C. 4> ‘ ) 

TYPE  1123 

1123  FCRMATC  NEUTRON  37  GROUP  MIN  IS  1.CE-5  EV;  PHOTON  21  CROUP  MAX’ 

1  'IS  1.0EM  EV‘) 

ACCEPT  1C11.DLBWI 

C 

C  ACCEPT  BOUNDS  ON  NORMALIZING  INTEGRAL 

C 

IF( INPUT2  .EQ.  C)CO  TO  212 
TYPE  1  Cl  3 

1  Cl  3  FORMATC  WHAT  IS  UPPER  BOUND  CN  NORMALIZING  INTEGRAL  IN  EV:’ 

1  '  <E10.4> ’ ) 

TYPE  1122 
ACCEPT  1011, DUENI 
TYPE  1014 

1  Cl 4  FCRMATC  WHAT  IS  LCWER  BCUND  ON  NORMALIZING  INTECRAL  IN  EV:  ' 

1  ' <E1C. 4>  '  ) 

TYPE  1123 
ACCEPT  1 C 1 1 , DLBNI 
212  CONTINUE 

O 

C  SET  UPPER  EOUND  ON  WEIGHTING  INTECRAL 

C 

IXCP-NRC¥*1 
XDUB'DUEWI 
DO  151  1*1,  NRCW1 
XE( I )  *ENGP( I ) 

151  CONTINUE 

CALL  UPBND 

IUBGWI-IXUBC 

UEFWI*XUBF 

IUBFGW-IXUBFG 

UBFEW-XUBFE 

C 

C  SET  UPPER  EOUND  CN  NORMALIZING  INTFGhAL 

C 


IF( INPUTS  .EC.  0 )GC  TC 

XDUB'DUBNI 

CALL  UPEND 

IUBFGN-ixUEFG 

UEFNI*XUBF 

IUBGNI-IXUBG 

UBFEN-XUEFE 


213 


21 


CONTINUE 


SET  LOWER  BOUND  WEICHTINC  INTEGRAL 


XDLB-DLBWI 

CALL  LOVBND 

ILBGWI-IXLBG 

LBFVI-XLBF 

ILBFCW-IXLBFG 

LBFEW-XLBFE 

LBFEW*XLEF~ 


SET  LOWER  BOUND  ON  NORMALIZING  INTEGRAL 


IF ( INPUT2  .EC.  C)GO  TO  214 

XCLB-CLBNI 

CALL  LCWENC 

IL3FCN-IXLBFC 

LBFNI-XLEF 

I LECNI *IXLEC 

LEFEN=XLBFE 

CONTINUE 


EVALUATION  CF  WEIGHTING  INTEGRAL 


DC  HC  I-1..NROW 
FRCTWti; =C . C 
FfiCTN(I',*C.C 
CONTINUE 


DETERMINE  DELTA  TERMS  FOR  INTEGRALS  I.L.  LE  OR  DU 
PLACE  ALL  LATA  IN  THE  FORM  CF  FKI'E.  IRFECARDLESS 

Ir(NSPT)  122,125,124 
DC  125  1=1 .NROW 
:SITA(I)*DE(I) 

FHIOT-  ,'I '  »FHI  T  ( I ) 

CONTINUE 
GO  TO  126 
DO  127  1*1 ,NRO¥ 

DELTA ( I )  *1 .0 

PHIENR (I)*PHITfI)/DE(I) 

CONTINUE 
GO  TO  126 
DO  126  I  *  1 , NRCW 
DELTA  f I }*IU ( I ) 

FHIENRf I ) *PHIT  v I  )*DU ( I  )/DE( I ; 

CONTINUE 


EVALUATION  OF  WEIGHT INC /FOLD INC  INTEGRAL 


SUM-O.C 

DC  1C,C  I  *IUEGWI ,  I  LEG  VI 


rr  »t  ■  r 


'.  ,*V  .  .  •.  *--*y 


FRCTW (I ) -PHIT(  I  )*DELTA  ( I  '*FNCT(  I 
SUM-SUM*FRCTW(I) 

150  CONTINUE 

UTAIL*PHIT( IUBFGW)*DELTA( IUEFCW )*PNCT( IUBFGW )*UBFWI 
LTAIL“PHIT(ILBFGN)*DEiTA( ILBFGW  J#FNCT( ILBFGW )*LBFWI 
IFdUBCWI  .BQ.  IUBFGW )G0  TO  230 
FRCTW ( IUBFGW ) -UTAIL 

23C  IFdLBCWI  .EQ.  ILBFCW)CO  TC  231 

FRCTWdDBFCW)*LTAIL 
231  WINT»SUN*UTAIL*LTAIL 


1 


L 

C  EVALUATION  OF  NORMALIZING  INTEGRAL 

C 

IF( INPUT 2  .EQ.  C )GC  TC  215 
SUM=O.C 

DO  1  11  I-IUECNI.ILBCNI 
FRCTN( I ) -PHIT( I )*EELTA ( I ! 

SUK=SUM*FRCTN( I ) 

111  CONTINUE 

UTAIL  =  PHIT ;  IUEFGN  )*DELTA  ( IUEFCN  ;*UEFNI 

ltail*phit(ilbfgn)*eelta(ilefgn;*lffni 

IF(IUBGNI  .EC.  IUEFGN )GC  TO  232 
FRCTN( IUBFCN ) 'UTAIL 

232  IF( ILBCNI  .EC.  ILEFGNdC  TC  233 
FRCTN (ILBFCN )  *LTAIL 

233  NINT-SUM»UTAIL»LTAIL 

C 

C  EVALUATE  EFFECTIVE  FUNCTION 

C 

FNCEFF=WINT/NINT 
215  CONTINUE 

C 

C  SET  OF  CE  FOR  DELTA  E  IN  THE  FCLIEI  FUNCTION  VAFIATKN  FCUATICN 

C 

IF( INPUT'  .EC.  C )GG  TO  216 
IC-IUBFCW-1 
LG  217,1*1 ,IC 
DEW (I ) *0 . C 

217  CONTINUE 
IC*(NR0W*1 )-ILBFGW 
DO  218  1-1  .  IC 
L*(NR0W1  )-I 

DEW ( L) *0.0 

218  CONTINUE 

DC  219  I  *  I UBFGW , I LBFGW 

dew(i;=de(i; 

215  CONTINUE 

IFdUBFCW  .EC.  IUECWI  )CC  TC  236 
DEW ( IUBFGW )=DEW ( IUBFCW )*UBFEW 
*36  IF(ILBFCW  .EC.  IL-’CWdGC  TC  237 

dew(ilbfcw;=dew(ilefcw)«lbfew 

237  CONTINUE 


COMPUTE  ESTIMATED  VARIATION  IN  THE  FOLDED  FUNCTION 
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.'  %’  *.■  *."  \  ‘  ,'7/  ^  •  . 


4 


a 


*1 


4 


.1 

■ 


-I 


c 


SUX-0.0  '-j 

DO  220  I-1.NR0W 

PWPHI(I )=(FNCT(I )*DSF( I )*PHIENR(I )*DEW(I ) )#*2 

PWKNT ( I ) - ( PHIESR ( I ) *DFF ( I ) *FNCT ( I ) *  DEW ( I ) ) **2 
PWSUM ( I ) -PWPH I ( I ) *PVTHT ( I ) 

SUM-PWSUMd)*SUM 
CONTINUE 
VFFNCT=SUK**0.5 


SET  UP  OF  DE  FOR  DELTA  E  IN  THE  NORMALIZING  INTEGRAL  PORTION 
OF  THE  EXPECTED  VARIATION  IN  THE  WEIGHTED  FUNCTION 


I F ( INPUT2  .EQ.  0 )C0  TO  216 

IC-IUBFCN-1 

DO  228  1  =  1  .  IC 

DEN ( I ) *C. 0 

CONTINUE 

IC-NROW-ILEFGN 

DO  22 1  1*1, IC 

L*(NR0W*1 )-I 

DEN  d)  *0.0 

CONTINUE 

DO  222  I*IUBFGN, ILEFGN 
DEN( I )*DE(I ) 

CONTINUE 

I F { IUBFCH  .EC.  IUBGNI)CO  TO  23 4 
DEN ( IUBFGN )*DEN ( IUBFCN )*UBFEN 
IF(ILBFCN  .EC.  ILBCNI)CO  TO  235 
DEN( ILBFCN)»DEN(ILBFCNI )*LBFEN 
CONTINUE 


CONFUTE  EXPECTED  VARIATION  IN  WEIGHTED  EVNCTICK 


SUM*0. 0 

DO  223  1=1 ,  NROW 

TERM1 (I5-FNCT(I)*DEW(I)/NINT 

TERM2 ( I ) -WINT*DEN ( I ) / (NINT**2 ) 

XXX  * ( TERM  1(1) -TERM2 ( I ) ) 

PWFPHI ( I ) * (XXX*PHIENR ( I )*DSF( I ) )**2 

TERM 3 ( I )*PHIENR(I )*DEW(I )/NINT 

PWFFNT (I ) - ( TERM3 ( I )*FNCT ( I } «DFF (I ) ) **2 

SUM-PWFPHI ( I ) ♦ PWFFNT ( I ) *SUM 

CONTINUE 

VWFNT-SUM**0.5 

CONTINUE 


OPEN  OUTPUT  FILE 


OPEN (UN IT *1 , NAME-FI LEA, TYPE* ‘NEW  1 ) 


OUTPUT  ALL  DATA 


WRITEO  ,1124) 


TYFE  1124 

1124  FORMATC  GROUP  ENERGY  PHI  DELETA  FUNCTION  ') 

DO  129  1*1 ,NROW 

WRITE ( 1,1125)1. ENGP( I ) , PHIT( I ) , DELTA ( I ) , FNCT ( I ) 

TYPE  1125,  I ,ENGP(I ) ,PHIT(I ),DELTA (I ) , FNCT( I ) 

1125  FORMAT ( 15, 4E11 .4) 

129  CONTINUE 
WRITE(1 ,1126) 

TYPE  1126 

1126  FORMATC  GROUP  ENERGY  FRACTION  NUMERATOR  FRACTION  DENOMINAT ' 

1  'OR') 

DO  130  1-1 ,NROW 

WRITE( 1 , 1 1 27 )  ICRCUP{ I ) , ENGP( I ) , FRCTW (I ) , FRCTN ( I ) 

TYPE  1127, IGROUP( I ),£NCP(I ) , FRCTW (I ) , FRCTN ( I ) 

1127  F0RMAT(I5 , El  1 ,4,8X,E1 1 .4, 12X.E1 1 . 4  ) 

130  CONTINUE 
TYPE  1128 
WRITE ( 1,1128) 

1128  FORMAT (2CX , '  NUMERATOR  INTEGRAL  DENOMINATOR  INTEGRAL'; 

WRITEC 1,1129)  IUBGWI , IUBGNI 
TYPE  1129,  IUBGWI, IUECNI 

1129  FORMATC  UPPER  BOUND  GROUP  ' ,  10X ,  13, 1 8X ,  13 ) 

WRITE(1 ,1130)  ILBGWI , ILBGNI 
TYPE  11 30, ILBGWI .ILBGNI 

1130  FORMATC  LOWER  BOUND  CROUP  ' ,  1CX,  13, 18X,  13) 

WRITE( 1 , 1131  )IUBFGW,IUBFCN 
TYPE  1131, IUBFGW, IUEFGN 

1131  FORMAT( '  FRACTION  GROUP  UE  ' , 1 0X , 13, 18X , 13 ) 

WRITE( 1 ,11 32)ILBFGW, ILBFGN 
TYPE  1132,  ILBFGW, ILBFGN 

1132  FORMATC  FRACTION  CROUP  LB  ' ,  1 CX ,  13 , 1 8X ,  13 ) 

WRITE(1,1133)UBFWI,UBFNI 
TYPE  1133,UBFWI,UBFNI 

1133  FORMATC  UPPER  BOUND  FRACTION', 9X. El  1  .3,1  IX, EH  .3) 

TYPE  1134,  LBFWI.LBFNI 
WRITEO, 1134)  LBFWI.LBFNI 

1134  FORMATC  LOWER  BOUND  FRACTION ' , 9X , El  1  . 3 , 1 1 X , El  1 . 3 ) 

TYPE  1135,  DUBWI , DUBNI 
WRITE ( 1,1135)  DUBWI, DUBNI 

1135  FORMATC  UPPER  BOUND  ' , 1 6X.E1 1 . 3 , 1 ’ X.E1 1 .3 ) 

TYPE  1136,  DLEWI , DLBNI 
WRITEO  ,11  36  )DLBWI ,  DLENI 

1136  FORMATC  LOWER  BOUND', 16X, Ell. 3. 11X, Ell. 3) 

TYPE  1137,  WINT.NINT 
WRITE ( 1,1137)  WINT.NINT 

1137  FORMATC  INTEGRAL  '  ,20X , El  1 . 3,  "  X ,  El  1 . 3 ) 

TYPE  1138,  FHCEFF 
WRITEO  ,1138)  FHCEFF 

1 1  38  FORMATC  WEIGHTED  FUNCTION  ' ,  1 7X ,  El  1 . 3 ) 

TYPE  1208.VFFNCT 
WRITE (1 ,1208),VFFNCT 

1208  FORMATC  VARIATION  ON  FOLDED  FUNCTION  \9X,E11.3) 

TYPE  1209,  VWFNT 
WRITEO, 1209), VWFNT 

1209  FORMATC  VARIATION  OF  WEICHTED  FUNCTION  \7X.E11.3) 

TYPE  121C 
WRITEO  ,1210) 

1210  FORMATC  CROUP  DELTA  FNCT  DELTA  SPECT  FHIFNR ' ) 

DC  224  1*1 ,NROW 

TYPE  1211 , 1 GROUP ( I ) , DFF ( I ) , DSF( I ) , FHIENR ( I ) 

WR ITE ( 1 , 1 2 1 1 ) , IGROUPf I ) , DFF ( I ) , DSF ( I ) , PH IENR ( I ) 

1211  FCRMAT( 15 , El  1 .4 ,8X , El  1 .4 , 1 2X , El  1 . 4 ; 


CONTINUE 
TYPE  1212 
WRITE(1 .1212) 

FORMATC  GROUP  DEW  PWPHI  PWFNT‘1 

DO  225  I-1.NROW 

TYPE  1213, IGROUP ( I ) , DEW(I ) , PWPHI ( I ) , PWFNT ( I ) 

WRITE( 1 , 1 21 3 ) , ICROUP( I ) , DEW ( I ) , PWPHI ( I ) , PWPNT( I ) 

FORMAT(I5,E1 1 .4, OX. El  1 .4, 12X, El  1 .4 1 

CONTINUE 

TYPE  1214 

WRITE(1 ,1214) 

FORMATC  GROUF  TERM!  TERM2  TERM?') 

DO  226  1-1 , NROW 

TYPE  1215, ICROUP( I ) , TERM1 ( I ) ,T£RM2( I ) , TERM3 ( I ) 

WRITE( 1 , 1215). IGROUP (I), TERM1 ( I ) ,TERK2( I ) ,TERM3( I ) 

FORMAT( 15 , El  1 . 4 ,8X , El  1  -  4 , 1 2X , E1 1 .4 ) 

CONTINUE 
TYPE  1216 
WR I TE ( 1 ,1216) 

FORMATC  GROUP  DEN  PWFPHI  PWFFNT ’ ) 

DO  227  1-1 , NROW 

TYPE  1217,  IGROUPd  ) ,  DEN ( I ) ,  PWFPHI  ( I  )  ,  PWFFNT ( I  ) 

WRITEO  ,  121 7),  IGROUPd  ),DEN{  I),  PWFPHI  (I ) ,  PWFFNT  (  I ) 

FORMAT( 15 , El  1 .4,8X,E1 1 .4.12X.E1 1  .4) 

CONTINUE 
CLOSE(UNIT-1  ) 

CONTINUE 

STOP 

END 


SUBROUTINE  UPBND 


SUBROUTINE  UPBNL 


THIS  SUBROUTINE  CALCULATES  THE  UPPER  BOUND  ON  INTEGRALS  IN  WHICH 
THE  INPUT  DATA  IS  IN  CROUPS  BUT  THE  LIMIT  IS  BETWEEN  THE 
BOUNDS  OF  A  CROUP 


COMMON  ICP.DUB, E(38 ) , IUBG, IUBFC , UBF 
COMMON  DLB , ILBC , ILBFG , XLBF, UBFE, XLBFE 


IGP  -  NUMBER  OF  GROUPS 
DUB  -  DESIRED  UPPER  BOUND 

E(38)  -  ENERCY  GROUP  SPECTRUM  (UPPER  ENERCY  BOUNDS) 

IUBG  -  UPPER  BOUND  GROUP 

IUFBC  -  UPPER  BOUND  FRACTIONAL  GROUP 

UBF  -  UPPER  BOUND  FRACTION  OF  GROUP  FOR  PHI*DE  OR  PH  I • DU 

UBFE*  UPPER  BOUND  FRACTION  OF  ENERCY 

DLB  -  DESIRED  LOWER  BOUND 

ILBG  -  LOWER  BOUND  GROUP 

ILBFG  »  LOWER  BOUND  FRACTIONAL  CROUP 

XLBF  -  LOWER  BOUND  FRACTION  OF  GROUP  FOR  PHI*DK  DOR  PHI*DU 
XLBFE  -  LOWER  BOUND  FRACTION  OF  ENERGY 


LOOK  TO  SEE  IF  UPPER  BOUND  LIES  ON  A  BOUNDARY 


o  o  o  o  o  o 


1 


4 


ISET-0 

DO  1  1-1, ICP 
ISET-ISET*1 

IF (DUB  .EQ.  E(I))  CO  TO  100 
COHTIIUE 
CO  TO  101 
100  IUBG-ISET 
IUBFG-ISET 
UBF-0.0 
CO  TO  105 

C 

C  IF  UPPER  BOUND  LIES  BETWEEN  CROUPS,  FIND  FRACTION  OF  GROUP 

C 


1C1  ISET-0 

DC  103  1-1 ,ICP-1 
ISET-ISET-1 

IF(DUB  .CT.  E(I ) )  GO  TO  104 
1C3  CONTINUE 

104  IUBC-ISET 
IUBFC-ISET-1 

DELTAE*E( ISET-1 )-E(ISET) 

NOTE:  E( ISET-1)  IS  THE  HIGHER  ENERGY 

PARTE-DUB-E(ISET) 

XI -E(ISET-1 ) 

X2*E(ISET) 

XDELTA-ALOG ( XI )-AL0G(X2) 

XPART- ALOC ( DUB ) -ALOG ( X2 ) 

UBF-SQRT ( ( XPART/XDELTA ) **2 ) 

UBFE-PARTE/DELTAE 

NOTE:  IF  DUB-E(ISET)  THEN  DUB-E( ISET ) *0  AND  UBF-0 

105  TYPE  1 0C2 

1CC2  FORMAT ( '  TEMP  AT  SUBROUTINE  FORMAT  1002') 

TYPE  1000 

1000  FORMATC  IUBC,  IUBFC.UBF.DELTAE, PARTE' : 

TYPE  1001 , I UBG , IUBFG , UBF , DELTAE , PARTE 

1001  F0RMAT(2I5,3E12.4) 

RETURN 

END 

C 

C  SUBROUTINE  LOWBND 

C 

SUBROUTINE  LOWBNC 
C 

C  THIS  SUBROUTINE  IS  FOR  CALCULATING  THE  LOWER  BOUND  ON  INTEGRALS 

C  WHERE  THE  INPUT  DATA  IS  IN  GROUPS  AND  THE  CHOSEN  BOUND  IS  BETWEEN 

C  THE  BOUNDS  OF  A  CROUP 

C 


COMMON  ICP, CUB, E(38),IUEG, IUBFG, UBF 
COMMON  CLB , I LBG , I LBFC , XLBF , UBFE , XLPFE 
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t  ICP  -  NUMBER  OF  GROUPS 

C  DUB  -  DESIRED  UPPER  BOUND 

C  £(38)  -  ENERGY  GROUP  SPECTRUM  (UPPER  ENERCY  BOUNDS) 

C  IUBC  -  UPPER  BOUND  GROUP 

C  IUFBG  -  UPPER  BOUND  FRACTIONAL  CROUP 

C  UBF  -  UPPER  BOUND  FRACTION  OF  CROUP  FOR  PHI*rE  CR  PH  1*1:1! 

C  UBFE-UPPER  BOUND  FRACTION  OF  ENERCY 

C  DLB  •  DESIRED  LOWER  BOUND 

C  ILBC  *  LOWER  BOUND  CROUF 

C  ILBFC  -  LOWER  BOUND  FRACTIONAL  CROUP 

C  XLBF  ■  LOWER  BOUND  FRACTION  OF  CROUP  FOR  PHI*DE  OR  PHI*DU 

C  XLBFE-LOWER  BOUND  FRACTION  OF  ENERCY 

C 

c 

c  CHECK  TO  SEE  IF  THE  LOWER  BOUNDRY  LIES  ON  A  BOUNDRY 


ISET-C 

DO  1  1-1 , ICP* 1 
ISET-ISETM 

IF(DLB  .EQ.  E ( I ) ) GO  TO  IOC 
1  CONTINUE 

GO  TO  101 

IOC  ILBC-ISET-1 

ILBFC-ISET-1 
XLBF-0.0 
GO  TO  105 


l 

c  IF  LOWER  BOUND  LIES  IN  A  BOUNDRY  OF  A  CROUP  FIND  THE  FRACTION  CP 

C  OF  THAT  CROUP 

C 

101  ISET-0 

DO  103  1-1 ,  ICP 

IF(DLB  .CT.  E(I))G0  TO  104 

ISET-ISET-1 

103  CONTINUE 

104  ILEG-ISET-1 
ILBFG-ISET 

DELTAE«E( ISET-1 )-E( ISET ) 

PARTE-E( ISET-1 )-dlb 
X 1 -E( ISET-1 ) 

X2-E(ISET) 

XDELTA-ALOG (XI )-AL0C(X2) 

XPART*AL0C(X1 )-ALOG ( DLB ) 

XLBF -SORT ( (XPART/XDELTA )**2 ) 

XLBFE-PARTE/DELTAE 

105  TYPE  1002 

1002  FORMATC  TEMP  FROM  SUBROUTINE  LOWBND  FORMAT  1002’) 

TYPE  1000 

1000  FORMATC  ILBG ,  ILBFG , XLBF, DELTAE,  PARTE ' ) 

TYPE  1001,  ILBC, ILBFG, XLBF, DELTAE, PARTE 

1001  FORMAT (21 5 , 3E1 2 . 4 ) 

RETURN 

END 
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